Dispersion interactions between semiconducting wires 
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The dispersion energy between extended molecular chains (or equivalently infinite wires) with non-zero band 
gaps is generally assumed to be expressible as a pair-wise sum of atom-atom terms which decay as R^^. Using 
a model system of two parallel wires with a variable band gap, we show that this is not the case. The dispersion 
interaction scales as z^^ for large interwire separations z, as expected for an insulator, but as the band gap 
decreases the interaction is greatly enhanced; while at shorter (but non-overlapping) separations it approaches a 
power-law scaling given by z^', i.e. the dispersion interaction expected between metallic wires. We demonstrate 
that these effects can be understood from the increasing length scale of the plasmon modes (charge fluctuations), 
and their increasing contribution to the molecular dipole polarizability and the dispersion interaction, as the 
band gaps are reduced. This result calls into question methods which invoke locality assumptions in deriving 
dispersion interactions between extended small-gap systems. 

PACS numbers: 34.20.Gj 73.22.-f 31.15.ap 



INTRODUCTION 



the atom-atom R interaction. 
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Conventionally, the dispersion interaction between two 
systems is formulated in terms of correlated fluctuations 
described by local (frequency-dependent) polarizabilities 
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which gives rise to the familiar atom-atom -C'^ R^i, 
teraction at leading order [ill • This form of the interaction 
has been used with a good measure of success in studies 
of many systems, from gases to solids, including com- 
plexes of biological molecules and organic molecules. 
These are typically insulators; that is, they have large 
HOMO-LUMO gaps. 

It has been known for a very long time that the additive 
atom-atom form of the dispersion does not hold very well 
for metallic systems. The classic case of an atom interact- 
ing with a thin metallic surface shows deviations from the 
additive law: the correct treatment of the metal results in 
a R^^ power law yl], but summing over atom-atom terms 
leads to a R'^ interaction. More recently, strong devia- 
tions from additivity have been demonstrated in the inter- 
actions of extended systems of metals or zero band-gap 
materials with at least one nano-scale dimension yHTfl. 

The semi-classical picture is useful for understanding 
the source of the differences in the metallic and insulating 
cases. In an insulator, electronic perturbations decay ex- 
ponentially with distance. We can therefore treat electron 
correlations as being local. At lowest order, these local 
fluctuations give rise to instantaneous local dipoles, and 
the correlation between these local dipoles gives rise to 
the atom-atom R^^ interaction. However, in a zero band- 
gap material, electronic fluctuations are long-ranged, par- 
ticularly if one or two of the dimensions are nano-scale 
IJSl]. It is these long -ranged fluctuations that give rise to 
the non-additivity of the dispersion and deviations from 



While the insulating and metallic cases are now well 
understood, little is known of the intermediate, semi- 
conducting case. The nature of the dispersion interac- 
tion between extended molecules with finite but small 
HOMO-LUMO gaps less than about 0.2 a.u. (5 eV) re- 
mains an open question. A large number of important 
nano-molecules fall into this category, such as carbon 
nano-tubes and the 'lander' -type molecules that are used 
as organic conductors. The electronic structure of mate- 
rials made of these molecules depends strongly on struc- 
tures they assume in the bulk, often via self-assembly. 
Consequently it is very important to understand exactly 
how these molecules interact. The dispersion interaction 
is the dominant source of attraction between these n- 
conjugated systems. For want of a clearer understand- 
ing, many studies assume the usual insulating case for 
the atom-atom R^^ form of this interaction. As we shall 
show in this article, this is both qualitatively and quanti- 
tatively incorrect for semiconducting molecules. 



A word of clarification; we use the term 'non-additive' 
here to describe deviations from the additive atom- 
atom picture of the dispersion interaction between two 
molecules. It is also commonly used to refer to the devi- 
ation from pair additivity seen in the interactions of three 
or more distinct molecules, but we are not concerned with 
such effects here, except for a brief comment in the Dis- 
cussion. 



INTERACTING WIRES USING HUCKEL 
(TIGHT-BINDING) THEORY 

The dispersion energy appears at second-order in inter- 
molecular perturbation theory and is formally expressed 
in terms of the exact eigenstates and eigenenergies of the 
non-interacting systems yj]. In a mean-field theory the 
dispersion energy between two subsystems (A and B), 
can be expressed as a sum over the single-electron wave- 
functions localised to each subsystem: 
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where i,j (a,b) are occupied (virtual) single-particle 
wavefunctions in either subsystem A or B, and e, is the 
eigenvalue of the /-th wavefunction. Assuming two par- 
allel wires, aUgned parallel to the x axis and separated 
by a distance z, the integral which appears above has the 
form: 
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and is the Coulomb interaction between the charge den- 
sity 4'*4'a due to excitation / ^ a in subsystem A with 
ip*jipb in B. 

In periodic systems, the wavefunction can be ex- 
pressed in Bloch form, i.e. i/'/x) = e'''''''uj{x), and so 
the co-density \l/*{x)ijj„(x) - e'^'^°~'^'^'^u*{x)Ua{x) can have 
a long-wavelength modulation, whose electrostatic field 
will not be well described by a multipole expansion at 
separations comparable to this length scale. Furthermore, 
in small-gap systems, these excitations have plasmon 
character and contribute significantly to E^^^ , as we dis- 
cuss below. 

Hiickel (tight-binding) theory gives us a convenient 
formalism for evaluating the above expression for inter- 
actions between two one-dimensional wires. We consid- 
ered a two-band model Hamiltonian of the form 



H - y (J3a^2ia2i-\ + /3' a' li+ian + h.c). 



(3) 



where yS,y6' are alternating bond-strengths between adja- 
cent sites, as would be encountered in a chain of (H2),, 
or a TT-conjugated polyene. We computed the interactions 
between two such parallel wires, each consisting of 2« 
identical atoms, equally spaced at intervals ofd, such that 
the unit cell is of length 2d and contains two atoms. As- 
suming periodic boundary conditions over a crystal cell 
of length 2dn, then there are two bands per wavevector. 
The single-particle wavefunctions and energies of sych a 
system can be found analytically. The band structure is 
given by 



and the set of wavevectors by 
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P - P' corresponds to a uniform wire with energy eigen- 
values given by e{k) - +2yS| cos(W)|, which in the limit 
of large « has a vanishing band-gap at half-filling (i.e. is a 
metal). The opposite limit (/?' = 0) corresponds to a chain 
of isolated dimers with energy eigenvalues e{k) - +/3, in- 
dependent of wavevector, and corresponds to the perfect 
insulator By varying the ratio /3' //3 between and 1, we 
can very conveniently probe the dispersion interaction in 
the intermediate (semiconducting) regime, with the band- 
gap given by AEg - 2(J3-/3'). The required integrals can 
be evaluated jsll as: 

{ij\r-,^\ab} = ^^/:o(Gz)F,,(G)F,i(^,+^,-^„-^i-G) 

(6) 
where [G] are reciprocal lattice vectors of the primi- 
tive unit cell and Yia{G) is the (analytic) Fourier trans- 
form of u*Ua obtained by assuming highly-localised ba- 
sis functions on each atom. Kq is the zero-th order 
modified Bessel function of the second kind and is the 
Fourier transform of the potential generated by a one- 
dimensional lattice at a field point at z from the lattice [9|- 
1111 . K() decays exponentially with increasingly large ar- 
guments and so typically only the 10 smallest reciprocal 
lattice vectors need to be included in the summation in 
order to obtain converged results. 

The calculations presented here use 8401 Monkhorst- 
Pack A:-points to sample the Brillouin zone, which cor- 
responds in real-space to wires consisting of 16802 sites 
per crystal cell. Such fine k-point sampling was required 
to obtained converged results with respect to system size. 

As shown in fig. [U the dispersion interaction between 
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two uniform wires ifi' Ifi - 1) varies as 
the whole range of z, in good agreement with previously 
obtained results via RPaII] and QMClltl for metallic 
wires. The interaction between chains of isolated dimers 
iP' - 0) varies as z^^, as would be expected from the con- 
ventional picture of dispersion interactions. In the semi- 
conducting regime, there is no unique exponent which 
characterises the dispersion interaction over all z- How- 
ever, at small separations z, it tends to the metallic be- 
haviour (z^^), whereas large z it tends to insulating be- 
haviour z^^ ■ 



INTERACTING Hj CHAINS USING SAPT(DFT) 

The principal drawback of the Hiickel model is the 
lack of electron correlation, and hence screening, within 
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FIG. 1: The second-order dispersion energy of parallel wires 
modelled using Huckel theory obtained from eq. m. The bond- 
ing interaction is of alternating strength and is controlled by the 
ratio /3'//3- The lines are numerical power-law fits to the small-z 
and large-z data points. In the metallic (p = P) and perfect insu- 
lating if}' = 0) cases, a single power-law fits the whole range of 
z, with the expected exponents -2 and -5 respectively. The in- 
termediate cases show a cross-over between these two limiting 
regimes as z is varied. 
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FIG. 2: Distorted (H2)„ -chains. The distortion parameter is de- 
fined as 77 = y/x. In all our chains we have chosen x = 1.4487 
a.u. 



ratio of the alternate bond lengths. The HOMO-LUMO 
gap of a hydrogen chain can be modified by simply intro- 
ducing a difference in alternate bond lengths. The distor- 
tion parameter 77 gives us a convenient way to control the 
electronic structure of the chain and allows us to study 
the interactions between insulating, semiconducting and 
(near) metallic chains in one framework. 



The frequency-dependent density response functions 
that appear in eq. ^ were evaluated with coupled Kohn- 
Sham perturbation theory, using the PBEO functional 
with a hybrid kernel consisting of 75% adiabatic LDA 
and 25% coupled Hartree-Fock. The accuracy of this ap- 
proach for the dispersion energy of small molecules sur- 
passes that of M0ller-Plesset perturbation theory and ri- 
vals that of coupled-cluster methods il3Lll4ll6lll7ll . Fur- 
thermore, there is no qualitative change in our results 
when the fraction of HF exchange is increased, so the 
results presented here are unlikely to be artifacts of the 
shortcomings of coupled Kohn-Sham perturbation the- 
ory ilSl 191 . All calculations used the cc-pVDZ basis 
which will result in a significant underestimation of the 
strength of the contribution to the dispersion energy that 
arises from transverse polarization, but should better de- 
scribe the contribution arising from the longitudinal po- 
larization. Since it is the longitudinal polarization that is 
most important in these systems, we expect our results to 
be qualitatively correct. 



each subsystem. For more realistic calculations we can 
use the symmetry-adapted perturbation theory based on 
density-functional theory 112^411 (SAPT(DFT)), where 
the dispersion energy is evaluated not via a sum-over- 
states but through a coupled Kohn-Sham formulation 
based on the density response functions of the interact- 
ing molecules 11 1511 : 
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where aAjsi^, r'; w) are the frequency-dependent den- 
sity response functions of the molecules which describe 
the propagation of a frequency-dependent perturbation 
in a molecule, from one point to another, within linear- 
response theory. The second-order dispersion energy is 
exact in this formulation if exchange effects are ne- 
glected. 

We have studied the interactions between two parallel 
finite (H2),, (fig.EJ chains with n - 16 and 32 and distor- 
tion parameters i] - 2.0, 1.5, 1.25 and 1.0, where rj is the 



Dispersion energies are displayed in fig. [3] and the ef- 
fective power laws for the physically important separa- 
tions between 6 and 20 a.u. are shown in table U together 
with HOMO-LUMO gaps. For insulating chains with 
an additive dispersion interaction we would expect two 
regimes determined by chain length L: for z « L we 
should expect the infinite chain result, i.e., the effective 
dispersion interaction should decay as z"^, and for z » L 
we should recover the usual z^^ power law. This is what 
we see with the chains with the largest distortion param- 
eter, 77 = 2.0, which exhibit large HOMO-LUMO gaps. 
As the distortion parameter 77 approaches unity and the 
HOMO-LUMO gap consequently decreases, the devia- 
tion from the additive insulating case becomes increas- 
ingly apparent, and for the longest of the chains consid- 
ered here the dispersion interaction is significantly en- 
hanced. For example at 40 a.u. (roughly half the chain 
length) an 77 = 1 chain has a dispersion interaction two 
orders of magnitude larger than an 77 = 2 chain. Addi- 
tionally, we see increasing finite-size effects: the power 
laws for the n = 16 and 32 chains, which were very sim- 
ilar for 77 - 2.0, are considerably different for 77 = 1.0. 
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FIG. 3: The second-order dispersion energy of pairs of parallel 
(H2),, chains. Energies have been normalized by the number of 
H2 units n. The solid lines are dispersion energies for chains 
with « = 32 and the dashed lines for chains with « = 16 (only 77 
= 1.0 and 2.0). 
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Power-law behaviour of £, fitted to the form z 
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TABLE I: 

in the region from 6 to 20 a.u. Beyond 600 a.u. all chains ex- 
hibit dispersion energies with the 7-* power law because of their 
finite length. The HOMO-LUMO gaps, t^Eg, have been calcu- 
lated from the Kohn-Sham eigenvalues. Energies are in atomic 
units. 
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0.370 


4.89 


0.366 


4.84 


1.5 


0.280 


4.58 


0.270 


4.50 


1.25 


0.202 


4.20 


0.183 


4.09 


1.0 


0.099 


3.52 


0.057 


3.17 



MULTIPOLE EXPANSION 

We can understand these effects in terms of the mul- 
tipole expansion. The multipole form of the dispersion 
energy is usually formulated in terms of (local) atomic 
polarizabilities. This leads to the usual R^^ atom-atom 
interaction (at leading order) and would not account for 
the anomalous dispersion power laws that we have ob- 
served in the chains. However the complete distributed- 
polarizability description is non-local: that is, it describes 
the change in multipole moments at one atom in response 
to a change in electrostatic fields at another. 

We obtain the multipole form of eq. O by expand- 
ing the Coulomb terms using the distributed form of the 
multipole expansion |r - r'|"' - Q"Tf^Qu where a and 
h denote (atomic) sites and t and u are multipole indices 
— 00 for the charge, 10, lie and \\s for the components 
of the dipole, and so on. Q] is the multipole moment op- 
erator for moment t of site a and r°f are the interaction 
tensors that contain the distance and angular dependence 
(see RefyJ] for details). Inserting this expansion in eq. 



(|7|l we obtain the multipole form for the dispersion en- 
ergy: 
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Here r,",* is the interaction function between multipole 
t on site a in subsystem A and multipole u on site b in 
B, while a""' are the frequency-dependent non-local po- 
larizabilities for sites a and a', and may be expressed in 
terms of the frequency-dependent density susceptibility 
asQ 

r f Q';(r)a{r,r'\tj)Q^:ir')dhdh'. (9) 
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There are no assumptions made in deriving eq. dHJ, other 
than that spheres enclosing the atomic charge densities 
on different molecules do not overlap. We have calculated 
distributed polarizabilities using the constrained density- 
fitting algorithm Ii20il . 

Eq. ([8]l involves a quadruple sum over sites and is 
therefore computationally demanding. To reduce this 
cost, we normally make a simplification by localizing 
the non-local polarizabilities. That is, the non-local po- 
larizabilities — the a"j°' {(jJ) with a + a' — are trans- 
formed onto one or other site using the multipole ex- 
pansion 112 ll I22II . This transformation is possible only if 
the non-local terms decay fast enough with inter-site dis- 
tance. If not, the multipole expansion used in the trans- 
formation diverges and the localization is no longer pos- 
sible. As we shall see, this is precisely what happens in 
the (H2),, chains as the distortion parameter r\ decreases. 

An important feature of the non-local polarizability de- 
scription is the presence of charge-flow polarizabilities — 
terms with t 01 u - 00 that describe the flow of charge in 
a molecule — which are usually small and die off quickly 
with distance. The lowest rank charge-flow polarizabil- 
ity is aJJooQ: if V" is the potential at site a, the change in 
charge at site a is given by AQl 
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y , a"" (V" -V") 
^fl "00,00'^ ^ ' '■ 

For a system with a large HOMO-LUMO gap the charge- 
flow terms are expected to be short-range. The charge- 
flow polarizabilities of the (H2)32 chain with 77 = 2.0 can 
be empirically modelled reasonably well with an expo- 
nential, that is. 
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where raa' is the intersite distance and y w 0.5 a.u. In this 
case, the charge-flow polarizabiUties drop by more than 
an order of magnitude within a few bonds. This is illus- 
trated in fig.m However, as we reduce the distortion pa- 
rameter 77 the charge-flow polarizabilities decay more and 
more slowly with site-site distance, until, at 77 = 1 .0, they 
span the entire length of the chain. Even for the 77 = 1 .5 
chain, these non-local charge-flow polarizabilities can no 
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FIG. 4: Matrix representation of tiie ciiarge-flow polarizability matrix ffg^ ^ for (H2)32 chains witii ;; = 2, 1 .5, 1 .25 and 1 (clockwise, 
from top left). Sites a and a' are represented along the x and y axes. Because these terms span a number of orders of magnitude of 
both signs we have plotted Indttg^Qol). Colour scheme: Black rectangles correspond to terms of order 1 a.u. and white rectangles to 
terms of order 10"^ a.u. 



longer be localized without incurring a significant error, 
but for the chains with r] - 1.25 and 1.0 localization re- 
sults in a qualitatively incorrect physical picture. 

Moreover the charge-flow polarizabilities contribute to 
the dipole-dipole polarizability in the direction of the 
chain, and this contribution increases dramatically as the 
band-gap decreases. For a 64-atom H2 chain with ?? = 2, 
the static dipole polarizability aw = o'liciic parallel to 
the chain, calculated as described above, is about 410 
a.u., of which 180 a.u. is contributed by charge-flows. 
When ;; = 1, so that the H atoms are equally spaced, 
a XX is about 11350 a.u., larger by about two orders of 
magnitude, and all of this increase is attributable to the 
charge-flow effects. 

We have calculated the dispersion energy using eq. 
([8]) with distributed polarizabilities, including terms to 
rank 1 . Higher ranking terms can be included, but these 
are not needed for the hydrogen chains. These results are 
presented in fig. |5] for the (H2)32 chains with 77 = 2.0 
and 1.0. First of all, consider the insulating chain (77 - 
2.0): the charge-flow polarizabilities alone severely un- 
derestimate the dispersion energy, but when terms up to 
rank 1 are included the agreement of eq. (O with the non- 
expanded SAPT(DFT) value for E^_^^ is excellent for all 
inter-chain separations shown in the figure. However, at 
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FIG. 5: The second-order dispersion energy calculated using 
the polarizabilities of (112)32 chains with distortion parameters 
r\= \ and 2. The contribution from the charge-flow polarizabil- 
ities, ffo^oo' '^ displayed separately from the sum of contribu- 
tions up to the dipole-dipole polarizabilities. 



77 = 1, the dispersion interaction is entirely dominated 
by the pure charge-flow (i.e. o'qJJ m,) terms, the higher or- 
der terms involving the dipole polarisabilities making a 
negligible contribution. Since the overall dispersion in- 
teraction is greatly enhanced as r\ is reduced, this implies 
that these pure charge-flow terms enhance the dispersion 



interaction in this limit. 

Why do the charge-flow terms resuh in the anomalous 
power laws? The lowest rank charge-flow terms, ffgooo' 
that appear in eq. (|8]l are associated with T-functions for 
the charge-charge interaction: T^q/Joo - ^ni' where Rat 
is the distance between site a on one wire and site b on 
the other If the wires are separated by distance vector 
R = (0, 0, z) and x^ and x^ are distance vectors for sites 
a and b along the chains then 'Rgb = R - (Xa - x^,). The 
dispersion energy arising from just the charge-flow terms 
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In contrast to the dipole-dipole polarizabilities, the 
charge-flow polarizabilities satisfy the sum-rule 
2fl' Q'J'ooC'^) - 0' which is a direct consequence of the 
charge conservation requirement: ^ aA(r,r'\w)d?r' - 0. 
This leads to cancellation between the charge-flow 
contributions to the total dispersion energy, and to an R^^ 
distance dependence at long range, but the cancellation 
is incomplete at short range, and terms in R", 2 < « < 5 
also occur To see how this arises, consider the following 
length scales: (1) L^ = I/T; which is a measure of 
the extent of the charge fluctuations determined by the 
exponential decay of the charge-flow polarizabilities 
assumed in eq. (fTOl i. and (2) L, the chain length. 

• z < Lc'. Here the charge-flow terms dominate 
and contribute R^^ terms to the dispersion energy. 
From fig. m we see that Lc is largest for the near- 
metallic wire for which we see the effects of these 
terms (fig.|5]l to z ~ 60 a.u. 

• Lc '& z < L: In this region the extent of the 
charge fluctuations is small compared with Rab and 
only those Ra'iy close to Rat are important. We can 
therefore expand Ra'f about Rai, using the multi- 
pole expansion. The leading term in this expansion 
is 

£f (00,00) =«--y—— 
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where we have used the charge-flow sum rule. So 
we see that the R^^ contributions sum to zero in 
this region. However, the higher-order contribu- 
tions are non-zero. 



• L «: z: In this limit both Rgj, and Ra'iy can be 
expanded in a multipole expansion about R = 



(0,0, zy.R-} = |R-(x«-Xi)r' = IR-x, 
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z" - jXgijZ^ 7 and likewise forTJ^,^,. Once again, 
the leading terms vanish because of the sum rule, 
leaving only the effective dipole-dipole contribu- 
tion: 
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This explains the large-z charge-flow contribution 
to the dispersion energy shown in fig.|5] Notice that 
for the near-metallic wires, this contribution to the 
total molecular C(, coefficient dominates that from 
the dipole-dipole polarizabilities, but for the large- 
gap wires the opposite is true. 

In short, the charge-flow polarizabilities give rise to the 
changes in power-law of the dispersion energy and also 
contribute to an enhancement of the effective Q coeffi- 
cient, applicable at long range. 



DISCUSSION 

The physical picture which emerges from both the 
Hiickel approach and the ab initio calculations is the fol- 
lowing. 

• In systems with a finite but small gap, spontaneous 
charge-fluctuations (plasmon modes in the infinite 
case) introduce a secondary length-scale interme- 
diate in size to the interatomic distance and the sys- 
tem size. This length scale grows as ffie gap gets 
smaller. 

• For separations z small compared with this length 
scale, the dispersion energy arising from the corre- 
lated fluctuations has metallic character 

• At z that are large compared with the length scale 
of ffie fluctuations, the dispersion can be described 
using London's dipole approximation, giving z"^ 
behaviour, but the magnitude of the fluctuations 
now depends strongly on the band gap, leading to 
orders of magnitude enhancement over the insula- 
tor case. 

• In small-gap systems these fluctuations give rise to 
a strong non-additivity in the polarizability. For ex- 
ample, the ratio of the longitudinal static polariz- 
abilities for the near- metallic and insulating (H2)32 



chains is 28. The dispersion energy is proportional 
to the square of the polarizabiHty, that is, 784. This 
is roughly the ratio of the dispersion energies for 
these two cases but only at separations z much 
greater than the chain length. Any attempt to ex- 
trapolate this result to shorter distances results in a 
severe overestimation of the dispersion energy. 

• This effect is not a consequence of retardation 
(we use the non-relativistic Hamiltonian) or damp- 
ing (charge-density overlap is negligible). It orig- 
inates from the complex behaviour of the non- 
local charge-flow polarizabilities. These are terms 
of rank zero that describe charge fluctuations in 
the system and are associated with a delocalized 
exchange-correlation hole l23l 12411 . This delocal- 
ization can be quantified using the localization ten- 
sor 11241 I25II and may give us a quantitative method 
for defining the charge fluctuation length-scale, L^. 
We are currently investigating this possibility. 

• We have demonstrated that both the change in 
power-law with distance and the enhancement of 
the dispersion energy can be understood using non- 
local polarizabiHty models containing charge-flow 
polarizabilities. 

• It should come as no surprise that these effects are 
also strongly dependent on system size (see fig|3]l. 

These results call into question theoretical methods 
that impose locality so as to scale linearly with system 
size (LCCSD(T), LMP2), or that approximate the dis- 
persion energy using a pair-wise -C(,R'^ interaction, as 
is done with dispersion-corrected DFT methods and em- 
pirical potentials. In the former, these effects can be in- 
cluded by extending the region of locality, though at the 
cost of losing linearity in scaling, but the latter methods 
should not be applied to systems such as these. Even DFT 
functionals with a non-local dispersion correction, such 
as the van der Waals functional of Dion et al. ll26ll and 
the more recent functional of Vydrov and Van Voorhis 
1 2711 are unlikely to contain the correct physics because 
of an implicit assumption of locality in the polarizabiHty. 
These functionals will include many-body non-additive 
effects between non-overlapping systems, but not the 
non-additive effects within each system such as those de- 
scribed here. On the other hand, methods based on the 
random phase approximation and quantum Monte Carlo 
should be able to describe the non-additive effects de- 
scribed in this paper if finite-size effects are kept under 
control. 

In systems containing carbon atoms, such as n- 
conjugated chains, the contributions from the core elec- 
trons will at least partially mask those from the more 
mobile n electrons, so it is possible that the changes in 



power-law of the dispersion interaction will not be as dra- 
matic as those we see in the (H2),, chains. But we should 
nevertheless expect a high degree of non-additivity aris- 
ing from the charge-flow terms. One indicator of this non- 
additivity is the non-linear dependence of the molecu- 
lar polarizabiHty on system size. This has been already 
demonstrated above and is further su ppo rted by the ex- 
perimental work of Compagnon et al. 112811 on the polariz- 
abilities of the fullerenes C70 and C60 which are found to 
be in the ratio 1.33 rather than the ratio 70/60 =1.17 that 
would be expected if the systems were additive. There- 
fore, if reliable Ce -models are to be constructed for ex- 
tended :7r-conjugated systems, these models will need to 
absorb the effects of the non-additivity of the charge- 
flow polarizabilities. That is, they should be tailored to 
suit the electronic structure of the system rather than be 
transferred from calculations on smaller systems. The 
Williams-Stone-Misquitta (WSM) procedure ll29l - l3lll is 
one such method that is capable of constructing effec- 
tive local polarizabiHty and dispersion models that ac- 
count for the electronic structure of the system. In fact, 
the non-local models presented in this paper are derived 
in the first step of the WSM procedure. We are currently 
investigating the behaviour of WSM dispersion models 
for a variety of carbon systems. 

Finally, strongly delocalized systems will also have 
important contributions to the dispersion energy from 
terms of third order in the interaction operator, that is, 
from the hyperpolarizabilities. We have not considered 
such terms in this paper. Furthermore, in ensembles of 
such systems there will be strong non-additive effects be- 
tween molecules. We are currently investigating the na- 
ture and importance of this type of non-additivity. 
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